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Abstract. For a given graph G with positive integral cost and delay on 
edges, distinct vertices s and t, cost bound C £ and delay bound 
D £ , the k bi-constraint path (fcBCP) problem is to compute k 
disjoint st-paths subject to C and D. This problem is known NP-hard, 
even when fc = 1 [3]. This paper first gives a simple approximation al- 
gorithm with factor-(2, 2), i.e. the algorithm computes a solution with 
delay and cost bounded by 2 * _D and 2 * C respectively. Later, a novel 
improved approximation algorithm with ratio (1 + /3, max{2, 1 + In ^'\) 
is developed by constructing interesting auxiliary graphs and employ- 
ing the cycle cancellation method. As a consequence, we can obtain a 
factor- (1.369, 2) approximation algorithm by setting 1 4- In ^ =2 and a 
factor-(1.567, 1.567) algorithm by setting 1 -I- /3 = 1 + In Besides, by 
setting /? = 0, an approximation algorithm with ratio (1, O(lnn)), i.e. an 
algorithm with only a single factor ratio O(lnn) on cost, can be immedi- 
ately obtained. To the best of our knowledge, this is the first non-trivial 
approximation algorithm for the fcBCP problem that strictly obeys the 
delay constraint. 
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1 Introduction 

In real networks, there are many applications that require quality of service and 
some degree of robustness simultaneously. Typically, the quality of service (QoS) 
related problem requires routing between the source node and the destination 
node to satisfy several constraints simultaneously, such as bandwidth, delay, 
cost and energy consumption. Nevertheless, in networks, some time-critical ap- 
plications also require routing to remain functioning while edge or vertex failure 
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occurs. A common solution is to compute k disjoint paths that satisfy the QoS 
constraints, and use one path as an active path whilst the other paths as backup 
paths. The routing traffic is carried on the active path, and switched to the 
disjoint backup paths while an edge or vertex failure occurs on the active path. 
However, for some time-critical applications even the time to discover failures of 
routing and restore data transmission in backup paths is too long for them. For 
such applications, packages are routed via k paths simultaneously, and the traffic 
is switched from failed paths to functioning paths if edge or vertex failures occur, 
such that routing can tolerate fc — 1 edge (vertex) failures. Therefore, given cost 
and delay as the QoS constraints, the disjoint QoS Path problem arises as below: 

Definition 1 For a graph G ~ {V, E) and a pair of distinct vertices s, t V , 
a cost function c : E , a delay function d : E , a cost bound C S Z^ 

and a delay bound D G Z^ , the fc-disjoint QoS Paths problem is to compute k 
disjoint st-paths Pi, . . . , Pk, such that ^ < C and d{Pi) < D for every 

i=l,...,k 

i = 1, . . . ,k. 

This problem is NP-hard even when all edges of G are with cost [5], which 
results in the difficulty to approximate the fc-disjoint QoS Paths problem. An 
alternative method is to compute k disjoint with total cost bounded by G and 
delay bounded by D ( equal to kD in Definition 1), and then route the packages 
via the paths according to their urgency priority, i.e., route urgent packages 
via paths of low delay whilst deferrable ones via paths of high delay of the k 
disjoint paths. Therefore, The disjoint bi-constraint path problem arises as in 
the following: 

Definition 2 (The k disjoint bi-constraint path problem, kBCP) For a graph 
G ~ (V, E) with a pair of distinct vertices s,t , a cost function c : E ^ , 
a delay function d : E ^ , a cost bound G G Z^ and a delay bound D G 
the fc-disjoint bi-constraint path problem is to calculate k disjoint st-paths 
Pi,..., Pk, such that J2 c{Pi) < G and J2 d{Pi) < D. 

i—l,...^k i— 

This paper will focus on bifactor approximation algorithms for the fcBCP prob- 
lem, which are introduced as below: 

Definition 3 An algorithm A is a bifactor (a, /3)- approximation for the kBCP 
problem, if and only if for every instance of kBCP, A computes k disjoint st- 
paths of which the delay sum and the cost sum are bounded by a* D and (3 * G 
respectively. 

Since a /^-approximation with the single factor ratio on cost is identical to a 
bifactor (1, /3)-approximation, we use them interchangeably in the text. 

1.1 Related work 

This fcBCP problem is NP-hard even when fc = 1 [1] . To the best of our knowl- 
edge, this paper is the first one that presents non-trivial approximation algo- 



rithms for the fcBCP problem formally. However, a number of papers have ad- 
dressed problems closely related to /cBCP, in particular the k restricted shortest 
path problem (fcRSP), which is to calculate k disjoint st-paths of minimum 

cost-sum under the delay constraint ^ d{Pi) < D. An algorithm with bifac- 

1=1,. ..,fc 

tor approximation ratio (2, 2) has been developed in [3] for general fc, while 
no approximation solution that strictly obeys the delay (or cost) constraint 
is known even when fc = 2. For a positive real number r, bifactor ratio of 
(1 + i, r(l + 2(i££i±i))(i + g)) and (1 + i, r(l + 2(i££I±i) )) have been achieved 
respectively in |10I3| for the case k = 2 and under the assumption that the delay 
of each path in the optimal solution of fcRSP is bounded by 

Special cases of this problem have been studied. When the delay constraint is 
removed, this problem is reduced to the min-sum problem, which is to calculate 
k disjoint paths with the total cost minimized. This problem is known polyno- 
mially solvable [TT]. Moreover, when fc = 1, the problem reduces to the single 
bi-constraint path (BCP) problem, which is known as the basic QoS routing 
problem [J and admits full polynomial time approximation scheme (FPTAS) 
[4l9j . Recently, the single BCP problem is still attracting considerable interests 
of the researchers. The strongest result known is a (1 -I- e)-approximation due to 
Xue et al [m. 

Additionally, when the cost constraint is removed, the disjoint QoS problem 
reduces to the length bounded disjoint path problem of finding two disjoint 
paths with the length of each path constrained by a given bound. This problem 
is a variant of the min-Max problem of finding two disjoint paths with the 
length of the longer path minimized . Both of the two problems are known 
to be A'P-complete [H], and with the best possible approximation ratio of 2 in 
digraphs [S], which can be achieved by applying the algorithm for the min-sum 
problem in |11I12| . Contrastingly, the min-min problem of finding two paths with 
the length of the shorter path minimized is NP-complete and doesn't admit K 
approximation for any K > 1 [511312) . The problem remains NP-complete and 
admits no polynomial time approximation scheme in planar digraphs [7] . 

1.2 Our techniques and results 

The main result of this paper is a factor-(l -t- /3, max{2, 1 -|-ln i}) approximation 
algorithm for any < /3 < 1 for the fcBCP problem. The main idea of the algo- 
rithm is firstly to compute fc-disjoint paths with delay-sum bounded by aD and 
cost-sum bounded by (2 — a)*C, where < a < 2 is a real number, and secondly 
to improve the computed fc paths by novelly combining cycle cancellation |10| 
and cost-bounded auxiliary graph construction [T3]. The key technique to prove 
the algorithm's approximation ratio is using definite integral to compute a close 
form for the sum of the cost increment during the improving phase. 

As a consequence of the main result, we can obtain a factor- (1.369, 2) approx- 
imation algorithm by setting 1 -I- In =2, and a factor- (1.567, 1.567) algorithm 
by setting 1-1-/3 = 1-1- In and slightly modifying our algorithm (to improve 
either cost or delay that is with worse ratio). Nevertheless, by slightly modifying 



Algorithm 1 A basic approximation algorithm for the fc-BCP problem 
Input: A graph G = (y,E), each edge e with cost c(e) and delay d(e), a given cost 
constraint C £ and delay constraint D £ Z^; 
Output: k disjoint paths Pi, P2 . . . , Pk- 

1. Set the new cost of edge e as 6(e) — + ^3^; 

2. Compute the k disjoint paths Pi, P2 . . . , Pfc in G by using Suurballe and Tarjan's 
algorithm [11112) . such that X^iLi X^eeP ^(^) '® minimized; 

3. Return Pi, P2 ...,Pk- 



our ratio proof, we show that an approximation algorithm with ratio ( 1 , O (In n) ) , 
i.e. an algorithm with single factor ratio of Oilnn) on cost, can be immediately 
obtained by setting /3 = 0. To the best of our knowledge, this is the first non- 
trivial approximation algorithm for the fcBCP problem that strictly obeys the 
delay constraint. 

We note that our algorithms are with pseudo-polynomial time complexity, 
since the auxiliary graph we construct is of size 0{C * n). However, by using the 
classic polynomial time approximation scheme design technique i.e. for any 

small e > setting the cost of every edge to in G before the construction 

of auxiliary graph, we can immediately obtain a polynomial time algorithm with 
ratio ((1 -I- ^) * (1 -I- e), max{2, 1 -I- In -g} * (H- e)). We shall omit the details due 
to the paper length limitation. 



2 An improved approximation algorithm for computing k 
disjoint bi-constraint paths 

This section will first present a simple approximation method for computing k- 
disjoint paths with delay-sum bounded by aD and cost-sum bounded by (2 — 
a) * C, where < a < 2 is a real number, and secondly improve the computed 
k paths by balancing the value of a and 2 — a. Though the presented simple 
algorithm is with worse ratio than that of the algorithm for A; = 2 in [TO] , it suits 
the improving phase better. 



2.1 A basic approximation algorithm 

Observing that the difficulty of computing fc-disjoint bi-constraint paths mainly 
comes from the two given constraints, the key idea of our algorithm is to deal 
with one new constraint B instead of the two given constraints C and D. Our 
algorithm firstly assigns a new mixed cost 6(e) = -I- to every edge in 
graph, and secondly computes k disjoint paths with the new cost sum bounded 
hy B = % + % =2. Note that the second step can be accomplished in polynomial 
time by employing the SPP algorithm due to Suurballe and Tarjan |11I12| . The 
detailed algorithm is as in Algorithm [TJ 



The time complexity and performance guarantee of Algorithm [T] is given by 
the following theorem: 



Theorem 4 Algorithm{^runs inO{km\ogi^ni n) time, and computes k- disjoint 
paths with delay-sum bounded by aD and cost-sum bounded by {2^a)*C , where 
< a < 2 is a real number. 

Proof. The main part of Algorithm [1] takes 0(fcm log _|_ m n) to compute k- 
disjoint paths by using Surrballe and Tarjan's algorithm |11I12| . and other parts 
of the algorithm take trivial time. Hence the time complexity of the algorithm 
is 0(fcm logi_|_m n). 

It remains to show the approximation ratio. To make the proof concise, we 
denote by OPT an optimal solution for the fc-disjoint BCP paths problem, and 
SOL the solution of Algorithm[TJ Obviously X^egOPT ^(^) — ^ holds. Then since 
the k disjoint paths is with minimum new cost, we have 



Assume the delay-sum of the algorithm is a times of d{OPT), then following Al- 
gorithm[T]0 < a < 2 holds. Therefore, we have J^egsol 

a + From Inequality a + < 2 holds. That is, c{SOL) < 

(2 - a)c{OPT) < (2 - a)C. This completes the proof. 

Note that a differs for different instances, i.e. Algorithm [1] may return a solution 
with cost 2 * c{OPT) and delay for some instances, while a solution with cost 
and delay 2 * d{OPT) for other instances. Hence, the bifactor approximation 
ratio for Algorithm [T] is actually (2,2). 

In real networks, the two given constraints may not be of equal importance, 
say, delay is far more important comparing to cost. In this case, applications 
require that the delay of the resulting solution is bounded by (1 + fi)D, where 
< /3 < 1 is a positive real number. Apparently, we could get an algorithm 
similar to Algorithm [1] excepting setting the new cost as 6(e) ~ + The 
ratio of the new algorithm is given as below: 

Corollary 5 By setting the new cost as 6(e) = (5^^^ + for a given real 
number < < 1 , Algorithm]^ returns k paths with delay-sum bounded by aD 
and cost-sum bounded by * C , where Q<a<l + (iisa real number. 

Therefore the ratio of the algorithm is(l + /3, 1 + ^). 

The proof of Corollary [S] is omitted here, since it is very similar to the proof 
of Theorem [TJ According to Corollary our algorithm can bound the delay- 
sum of the /c-disjoint path by {1 -\- P)D for any < /3 < 1, by relaxing the 
cost constraint to {1 -\- ^) * C. For example, if /3 = 0.01, then the bifactor 
approximation ratio of the algorithm is (1.01, 101). Thus, the algorithm decrease 
the delay of the /c-disjoint paths at a high price. In the next subsection, we shall 
develop an improved method that pays less to make delay-sum of the A:-disjoint 
paths bounded by (1 -I- (i)D. 




6(e) < 2. 



(1) 
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Algorithm 2 An improved algorithm based on cycle cancellation. 



Input: A graph G = (y,E), each edge e with cost c(e) and delay d(e), a given cost 
constraint C £ and delay constraint D £ , disjoint QoS paths Pi, P2 . . ■ , Pk 
computed by Algorithm [T] 

Output: Improved disjoint QoS paths Q\, Q2 ■ ■ ■ , Qk- 

1. If Eti'^(^'0 < (i + /3)^ ; 

then return Pi, P2 . . . ,Pk as Qi, Q2 ■ ■ ■ , Qk , terminate; 

2. Reverse direction of the edges of Pi, P2 . . . ,Pk in G , set their cost to a small 
positive real number < e < , and negative their delay; 

3. Compute cycle Oj with c{Oj) < C, d{Oj) < and ■^^-^ attaining minimum, by 
the method given in next section; 

/* Following clause 2 of Proposition [51 if ELi '^(-^0 > d{OPT) and ^^"=1 c(P) > 
0, there always exist cycle Oj with c{Oj) < C and d{Oj) < 0. */ 

4. Improve Pi , P2 . . . , P*, by adding the edges of Oj and removing the pairs of parallel 
edges in opposite direction; 

5. Go to Step 1. 



2.2 The improving phase 

To make the delay of the solution resulting from Algorithm [T] bounded by (1 + 
/3)D, our improving phase is, basically a greedy method, using the so-called cycle 
cancellation to improve the disjoint paths in iterations until a solution with the 
best possible ratio (1 + /3, max{2, 1 + In ^}) is obtained. The cycle cancellation 
method is an approach of using cycles to change the edges of the disjoint paths, 
which first appears in [10] and is derived from the following proposition that can 
be immediately obtained from fiow theory [T]: 

Proposition 6 Let Pi, P2 . . . , Pk and Qi, Q2 ■ ■ ■ , Qk be two sets of k disjoint 
st-paths in G, G he G excepting that all edges of Pi, P2 . . . , Pk are reversed, and 
O be a cycle in G. Then 

1. The edges of Pi, P2 . . . , Pk and O, excepting the pairs of parallel edges with 
opposite direction, compose k-disjoint paths; 

2. There exist a set of edge disjoint cycles Oi, . . . , Oh in G, such that the edges 
of Pi, P2 . . . , Pk and Oi, . . . , Oh, excepting the pairs of parallel edges with 
opposite direction, compose Qi, Q2 ■ ■ ■ , Qk- 

From the proposition above, it is obvious that there exists a set of cycles Oi, . . . , Oh 
that can improve k disjoint QoS paths Pi, P2 . . . , Pk to an optimal solution. 
However, it is hard to identify all the cycles Oi, . . . , Oh, so we employ a greedy 
approach to compute a set of cycles to obtain an approximation approach. The 
improving phase is composed by iterations, each of which computes a cycle and 
then uses it to improve Pi, P2 ■ ■ ■ , Pk- More precisely, to obtain a good ratio, 
the algorithm computes in iteration j a cycle Oj with minimized among 

the cycles in G. The layout of the algorithm is as given in Algorithm [51 



Following clause 1 of Proposition [51 Algorithm [2] will correctly return k dis- 
joint paths. It remains to show the cost and delay of the k disjoint paths is 
constrained as below: 



Theorem 7 The approximation ratio of Algorithm 2 is max{2, 1+ln ^})- 

Proof. For the case that X]i=i d{Pi) < (1 + P)D holds before the improving 
phase, the approximation ratio of Algorithm [2] is obviously (1 + /3, 2). 

It remains to show the ratio of the algorithm is (1 + /?, 1 + In i) for the case 

that J2i=i d{Pi) > (1 + I3)D. Assume that Algorithm [2] runs in h iterations, the 
key idea of the proof is to sum up the cost increment while using the cycle to 
improve the k disjoint paths in iterations, and show that the cost sum is bounded 
(by giving the cost sum a close form). 

Note that in the case, we have aD > J2i=i d{Pi) > (1 + P)D, so a > 1 + /3 
holds. Let AD = d{OPT) - d{SOL) > (1 - a)d{OPT) and AC = c{OPT) - 
c{SOL) < {a - l)c{OPT). Clearly, AD < and AC > hold. Let the cycle 
computed in the jth iteration be Oj, then since | attains minimum in Step 

3 of Algorithm H we have < ^^"^-r' ''''^'^ . That is, 



By summing up c{Oj) in h—1 iterations (excluding the last iteration), we have: 

d{0,) 



'J ) 

h-l h-1 



^.AD-Y^Udm 



Following the definition of Definite Integral, we have: 



> Mi = > --i d(Oi) < / -dx, 

U AD - ElZl dm fr{ AD - EiZi dm ' J AD 

(2) 

where the maximum is attained when d{Oj) = — 1 for every j. 

Algorithm [2] terminates when d{SOL) + J2i=i d{Oi) < (1 + li)D, so in the 
h - 1 iterations d{SOL) + YJIZI dm > (1 + I3)D holds. That is d{SOL) - 
D + Y!IZI dm > PD, and hence ~AD + Y!IZI d{Oi) >/?£»> holds. So we 
obtain a close form for the cost sum of the /i — 1 iterations: 



.4D-Et"i^<i(0,) 1 r-AD ^ .-AD ^ i^^i 

/ —dx = / —dx < / —dx = In = In — 

I AD X J-AD+Y.'^I^ d(Oi)) X JpD ^ jjD 

(3) 



At last, the cost increment in the hth iteration is bounded by c{OPT). So 
the final cost is c{SOL2) < (2 - a)C + Cln + C = C(3 - a + In where 
SOL2 is the solution resulting from Algorithmic] 

Let /(a) = 3-a + \n^. Remind that a < 2, so /'(a) = ^ - 1 > 0, /(a) 
is monotonous increasing on a, and attains maximum while a = 2. So we have 
c{SOL2) < (l + lni)c(OPr). 

Therefore, the cost of the output of Algorithm[2]is bounded by (1+ln ^)c{OPT), 
and delay bounded by (1 + f3)d{0PT). This completes the proof. 

From Theorem[71 by setting 1+ln = 2, we can immediately obtain an improved 
algorithm with best possible delay ratio under the same cost bound 2C. That 
is: 

Corollary 8 By setting 1 + In = 2, we have (3 ~ ^, and hence Algorithm\^ is 
now with a bifactor approximation ratio of (1 + i, 2) = (1.369, 2). 

For those applications in which delay and cost are of equal importance, by setting 
1 + In = 1 + /? and slightly modifying Algorithm [5] to improve either cost or 
delay that is of worse ratio, we can obtain an improved algorithm with ratio as 
in the following corollary: 

Corollary 9 // 1 + In ^ = 1 + /3, Algorithm\^ is with a bifactor approximation 
ratio of (1.567, 1.567). 

Now we consider the case that /3 = 0, i.e. the delay constraint is strictly satisfied. 

In this case, Inequality ^ in the proof of Theorem[7]will become ^j-^i ^ 

Corollary 10 When /3 = 0, Algorithm\^is with a ratio o/(l, O(lnn)). 

From CoroUarv llOl we can see that the price of obeying one constraint strictly 
is very high, i.e. it requires extra 0(ln n) times of cost. However, this is the first 
algorithm with logarithmic factor approximation ratio for the /c-BCP problem 
with strict delay constraint. 



3 Computing Cycle Oj with minimum 

Let G = iy, E) be G, excepting that the edges of Pi , P2 . . . ,Pk are with direc- 
tion reversed, cost sat to 0, and delay negatived. This section will show how to 
compute a cycle O with cost bounded by C and minimized in G. The key 
idea is firstly to construct an auxiliary graphs H{v) for each v where every cycle 
is with cost at most C, secondly to compute the cycle O' with minimum 
among all cycles in all H{v)s for each v € G, and thirdly to obtain cycle O with 
minimum ^Igl^jj in G according to O'. 



Algorithm 3 Construction of auxiliary graph H. 



Input: Graph G = {y,E), two distinct vertices s, t £ V , & cost c : e ^ Zq and a 
delay d : e on every edge e £ E, a cost constraint C and a delay constraint D; 

Output: Auxiliary graph H{v). 

1. For every vertex vi of V, add to H{v) vertices ,vf ; 

2. For every edge e — {vj,vi) G E, add to H{v) the edges 
wj", f;'^''^'^''"^ , . . . , ^i;^^'^''^' , , each of which is with cost c(e) and delay 

die); 

/*Note that d(e) can be negative in G = (V, E).*/ 

3. For alH = 2, . . . , C, add to H{v) backward edge (u', u^) with delay and cost 0, 
where a backward edge is an edge w-*) where i > j. 

/*H{v) contains backward edges, and hence cycles, only after adding the edges of 
Step 3.*/ 



3.1 Construction of auxiliary graph H{v) 

The algorithm of constructing the auxiliary graph H{v) is inspired by the method 
of computing a single path subject to multiple constraints [T3]. The full layout 
of the algorithm is as shown in Algorithm [3] (An example of such construction 
is as depicted in Figure [T|). 

Following Algorithm [31 every backward edge in the constructed auxiliary 
graph H{v) must contain vertex v^. Hence every cycle in H{v) contains at most 
one backward edge. On the other hand, following Algorithm [3] a cycle in H{v) 
contains at least one backward edge. Therefore, there exist exactly one backward 
edge in any cycle of H[v). Because H{v) \ {{v^ , v^) , . . . , (w*^, w^)} is an acyclic 
graph where any path is with cost at most C, we have: 

Lemma 11 Any cycle in H{v) is with cost at most C. 

Let 0{v) be a cycle in H{v), then following the construction of H{v), 0{v) 
apparently corresponds to a set of cycles in G. Conversely, every cycle containing 
V in G corresponds to a cycle in H{v). Based on the observation, the following 
lemma gives the key idea of computing a cycle O of G with minimized and 
cost bounded by G: 

Lemma 12 Let 0{vi) he a cycle with minimum f|Q|"'jj i'^ Hivi), and 0{v) he 
the cycle with minimum among the n cycles 0{vi), . . . , 0(w„). Assume O 

is a cycle with minimum in the set of cycles in G that correspond to 0{v). 
Then for any cycle O' in G with c{0') < G , < j holds. 

Proof. Suppose this lemma is not true, then there must exist in G a cycle, say O' , 
such that ^ > and c(O') < G hold. Then the cycle 0'{v) in H{v) that 

corresponds to O' is also with = ^f§7j < < ^(§(^; contradicting 

with the minimality of 0{v) in H{v). This completes the proof. 




(1)) 

edge with delay and cost e 

edge with delay and cost equal to its corresponding edge in G 

Fig. 1. Construction of auxiliary graph H{v = s) with cost constraint C = 6: (a) graph 
G; (b) auxiliary graph H{v = s). The cycle O = syts in G is exclude in the auxiliary 
graph H{s) as shown in (b), keeping the cost of k disjoint paths constrained by C = 6. 



3.2 Computing the cycle O with minimum -^^j^ 

The main idea of the algorithm to compute a cycle O with minimized in G 
is to compute the cycle O' with minimum ^[^tj among all cycles in all H{v)s for 
each w g G. Following Lemma [T^l the cycle O in G is the cycle with minimum 
among the cycles in G corresponding to all the computed O's. The detailed 
steps are as below: 

1 . For i = Iton 

(a) Construct H{vi) for 6 G by Algorithm [31 

(b) Compute cycle 0{vi) with minimum e(o("')) H(vi) by employing the 
minimum cost-to-time ratio cycle algorithm in [T]; 

(c) Select 0{v) with minimum jl^l^jj from the n computed cycles 0{vi), . . . , 0(w„); 

2. Select the cycle O with minimum among the cycles in G that correspond 
to 0{v). 

Clearly, the cycle O attains minimum in G. Besides, following Lemma [TT] 
we have c(0) < C. Therefore the cycle O is correctly the promised cycle. This 
completes the proof of the approximation ratio. 

4 Conclusion 

This paper gave a novel approximation algorithm with ratio (1 -|~ /3, max{2, 1 -|- 
In -^l) for the /cBCP problem based on improving a simple (a, 2— a)-approximation 
algorithm by constructing interesting auxiliary graphs and employing the cycle 
cancellation method. By setting (3 = 0, an approximation algorithm with bifac- 
tor ratio (1, O(lnn)), i.e. an 0(lnn)-approximation algorithm can be obtained 
immediately. To the best of our knowledge, it is the first non-trivial approxima- 
tion algorithm for this problem that obeys the delay constraint strictly. We are 
now investigating whether any constant factor approximation algorithm exists 
for computing a solution that strictly obey the delay constraint. 
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